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ABSTRACT 

Context. The presence of the moat flow around sunspots is intimately linked to the mere existence of sunspots. 

Aims. We characterize the moat flow (MF) and Evershed flow (EF) in sunspots to enhance our knowledge of sunspot structures and 

photospheric flow properties. 

Methods. We calibrated HMI synoptic Doppler maps and used them to analyze 3h time averages of 3 1 circular, stable, and fully 
developed sunspots at heliocentric angles of some 50°. Assuming axially symmetrical flow fields, we infer the azimuthally averaged 
horizontal velocity component of the MF and EF from 51 velocity maps. We studied the MF properties (velocity and extension) and 
elaborate on how these components depend on sunspot parameters (sunspot size and EF velocity). To explore the weekly and monthly 
evolution of MFs, we compare spots rotating from the eastern to western limbs and spots that reappear on the eastern limb. 
Results. Our calibration procedure of HMI Doppler maps yields reliable and consistent results. In 3h averages, we find the MF 
decreases on average from some 1000 + 200 m/s just outside the spot boundary to 500 m/s after an additional 4 Mm. The average 
MF extension lies at 9.2 + 5 Mm, where the velocity drops below some 180 m/s. Neither the MF velocity nor its extension depend 
significantly on the sunspot size or EF velocity. But, the EF velocity does show a tendency to be enhanced with sunspot size. On a 
time scale of a week and a month, we find decreasing MF extensions and a tendency for the MF velocity to increase for strongly 
decaying sunspots, whereas the changing EF velocity has no impact on the MF. 

Conclusions. On 3h averages, the EF velocity scales with the size of sunspots, while the MF properties show no significant correlation 
with the EF or with the sunspot size. This we interpret as a hint that the physical origins of EF and MF are distinct. 
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1. Introduction 

As complex magnetic formations, sunspots strongly affects the 
photospheric dynamics of the granular convective pattern of the 
plasma in which they are embedded. They suppress the upward 
propagation of heat in the magnetic core, seen as the cooler and 
therefore darker umbra, and they harbor the well-known penum- 
bral Evershed flow observed as a radially outwards-directed, 
horizontal flow in the photospheric layers. In photospheric ob- 
servations this leads to a blueshift and blueward asymmetry of 
the spectral line for penumbral filaments facing the center of 
the solar disk and a redshift and redward asymmetry for fila- 
ments facing the solar limb (e.g., |Evershed|19 09 ; Maltby 1964; 
|Schl ichenmai er et al.|2004| ). In general, the line shift, as well as 
the asymmetry, increases from the inner to the outer penumbra. 
In photospheric layers the Evershed flow is close to horizontal 
with maximum inclinations of 5° to 10°, exhibiting downward 



flows in the outer penumbra (e.g., Michard 1951 Balthasar et al. 



1996 



Schlichenmaier & Schmidt 2000; Franz & Schlichenmaier 



2009). Average horizont al velocities are in the range of 3-4 
km/s ( [Shine et"aL] |1994|), but can exceed 6 km/s or more on 



small scales (e.g., 



Roup pe van der Voort| 2002 ; 



Bellot Rubio 



et al.||2003a| ). The Evershed flow seems to stop abruptly at the 
white light boundary of the spot ( [Wiehr &" Degenhardt 1992; 
Schlichenmaier & Schmidt 1999]), which is in line with down- 
ward flows in the outer penumbra. Only a small fraction of the 
flow continues to flow into the magnetic canopy that surrounds 
the spot ( Reza eTet al.|2006| ). 

Whereas the Evershed flow (EF) is magnetized (e.g., 



Schlichenmaier & Collados 2002), the annular region around 



fully-fledged sunspots, called the moat cell, is largely nonmag 
netic. It harbors the moat flow (MF), a radially oriented outflow 
of plasma adjacent to the outer penumbral boundary. Embedded 
in these motions, small, moving magnetic features (MMFs) mi- 
grate away from the sunsp ot ( [Sheeley 1 1 969| 1 1 972} | Vrabecj 1971 



Harvey & Harvey] [1973| ). The existence of bo th the MF anc 
the EF depend on the presence of a penumbra (Vrabec p974 ). 
Consequently, the moat flow can only be observed fo r sunspot 
sides with a well-developed penumbra and lacks pores (Sobotka 
|etal.|1999| ). 

The moat flow develops immediately after the formation of 
a penumbra ( Pard on et al.|1 979 ; Schlichenma ier et al.|2010|) and 
was observed even after the decay of the penumbra ( [Verma et al. 
2012). The MF is present as a stationary flow through the spot's 
lifetime, which varies between several weeks to several months 
(see, e.g., |Solanki|2QQ3] ): [Sobotka & Roudier| ( [2QQ7l ) do not find 
any significant variations in the MF within half a day. During the 
initiation of the MF, the magnetic components in the vicinity of 
the sunspot are pushed to the periphery and leave the moat cell 
largely nonmagnetic with single MMFs. The extension of the 
moat field reaches between 10 Mm to 20 Mm from the penum- 
bral boundary for small sunspots and can be roughly twice the 
spot radius for larger spots (Brickhouse & Labonte||T988 ), but 



Sobotka & Roudier ( 2007J) find no such correspondence. 



The moat flow velocity ranges between 0.5 km/s to 1 km/s 
and can be seen by tracking the bright granular f eatures (e.g., 
Rimmele||1997| | Vargas Dommguez et al.|[2007t IBalthasar & 



Muglach|2010|) or Doppler shift measurements (Balthasar et al. 



|1996| ). That the MF consists of migrating granules, indicates that 
the MF is as nonmagnetic as granules, i.e., largely nonmagnetic. 
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Also helioseismic measurements have revealed the existence 
of the MF as an outflow extending up to 30 Mm with a maximum 
velocity of 1 km/s just next to the penumbral boundary ( |Gizon| 
et al.|2000] ) in the first 2 Mm of the solar surface. According to 



recent studi es ( [Sun et al.|19971|Gizon et al.|2009"l|Featherstone| 
et al.|20lT] ), the MF is also detectable in deeper layers, but has 



slo wer speeds compared to surface measurements. 

Cabrera Solana et al. (2006) suggest that a link exists 



between EF and MMFs. Magnetic velocity packages, called 
the Evershed clouds inside the penumbra ( [Shine et~aL 1994| 
Rimmele |l994 Cabrera S olana et al.|2007| ) propagate outwards 



to the extension of penumbral filaments and the moat region 
where they are embedded as MMFs. MMFs can travel from 
the penumbra into the vicinity of the sunspot ([S ainz Dal da & 



Martinez Pill et| [20051 [Sainz Dalda & Bellot Rubio||2008a|b| ). 
Also |Sch lichenmaier ( 2002j) proposes a scenario that is consis- 
tent with observed MMFs. In this way, a magneto-convective 
overshoot instability in an Evershed flux tube leads to the mi- 
grating feature in the MF regio n. There is still discordance about 
a link between the EF and MF. Vargas Dominguez et al. ( 2007| ) 
observed irregular sunspots and found that the moat flow is only 
present in radial extensions of penumbral filaments, but not per- 
pendicular to them. This indicates that the MF is an extension 
of the EF. However, this seems to be impossible since the EF is 
magnetized, while the MF is intrinsically unmagnetic, as men- 
tioned above. 

In this paper, we analyze the EF and MF and elaborate on a 
possible link between the two. To that end, we utilize Doppler 
shift measurements of HMI. In Sec. [2] we perform a thorough 
calibration of HMI Doppler maps. In Sec. [3] we describe the cri- 
teria for data selection and the method we used to analyze the 
flow properties. In Sec. [4] we present our results for the flow 
properties and elaborate on correlations between spot and MF 
properties and on how they change during the evolution of a 
spot. The statistics are based on 51 maps constructed from 31 
different spots. The data sets we used and the radial dependen- 
cies of the flow velocity of all spots, which form the basis of our 
analysis, are given in the Appendix. In Sec. [5] we discuss the 
present understanding of MFs in the context of our results and 
summarize our results and conclusions. 



2. Calibration of HMI Doppler maps 

Data: Our work is based on synoptic 720s intensity maps and 
Doppler map^]in the Fel spectral line at 6173.3 A from the 
JSOC webpage recorded by the Helioseismic Magnetic Imager 
(HMI) of the Solar Dynamics Observatory (SDO) up to maxi- 
mum values of ± 6.5 km/s and a spatial resolution of « 0.5"/px. 



Calibration by subtracting systematic components: To yield 
undisturbed flow velocities relative to the solar surface, we need 
to define a rest frame (v = 0). To this end, we construct a time 
averaged velocity map, (v(x,y)) t , in which the effects of granu- 
lation, oscillations, and supergranulation are removed. This av- 
erage is composed of three systematic, large-scale components: 



(v(x,y)) t = v c iv(r) + VrotGf, /, B ) + v res U, y) 



The first term is the center-to-limb variation in the convective 
blueshift, v c i v . It is radially symmetrical and is several hundred 
m/s, depending on the distance, r, (or heliocentric angle 6) to disk 



center. The second term is the differential rotation, v rot . It is ax- 
ially symmetrical with velocities up to ±2000 m/s with respect 
to the meridian. The values depend on the distance d from the 
meridian, the latitude / and the inclination Bo of the rotation axis, 
which varies annually between ±7.27°. The third term is a non- 
symmetrical residual, v res , with velocities up to ±150 m/s, which 
contains instrumental effects and other systematic flow fields of 
smaller magnitude. The large-scale meridional flow with veloci- 
ties of less than 20 m/s can be neglected for our sunspot studies, 
and if existing, it is included in these three components. In the 
next four sections we determine the four terms of Eq. [T] 

2.1. Construction of a time-averaged velocity map: (v(x,y)) t 

In velocity maps, supergranules produce a redshift on the limb 
side and a blueshift on the opposite side when they are off disk 
center. Because a supergranule migrates across the disk its aver- 
age velocity signal vanishes when it is next to the equator, since 
blueshifts are eliminated by the following redshifts. However, 
close to the poles, the blueshifts and redshifts are separated in 
latitude such that a single supergranule leaves noticeable traces 
in the Doppler map. To diminish these traces a long time se- 
ries is necessary. In our case we have averaged 720s Doppler 
maps of 61 entire days (a total of 7320 maps ) between June and 
December 2010, shown in the left hand panel of Fig.[T] 
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Fig. 1. — Left panel: Averaged velocity map, (v(x,y)) ? =6idays> 
with maximum blueshifts (eastern limb) and redshifts (western 
limb) with some 2 km/s. The axes are given in arcsec. — Right 
panel: Reduced velocity map, (v(x,y)) t - v c \ Y . 



In order not to spoil the time average by active regions, we 
have masked regions according to the darkening in the con- 
tinuum intensity and replaced the content of the mask in the 
Doppler map by the velocities on the opposite side of the he- 
liocentric equator. Thanks to the solar activity minimum, this 
method was appropriate. For consistency, we normalized the 
solar radius of all Doppler maps and eliminated the respective 
observer motioiQ relative to the Sun. Since we favored time 
periods with inclination Bo ~ 0°, the obtained velocity map 
(v(x, y)) t =6i days should serve as a good basis for modeling the 
systematic components named in Eq.[T] In (v(x,y)) t the velocity 
at disk center is about 50 m/s. 



(1) 2.2. Radially symmetrical velocity component: v dv (r) 



As seen in the left hand panel of Fig. [T] the LOS velocity of 
the solar rotation, which should be zero for the meridian and ax- 
ially symmetrical with respect to the centerlines, is deformed. 



1 A velocity map of the Sun according to spectral line shifts due to 
the Doppler effect. 



2 Keywords for observer motion in the FITS file header: OBS -VR 
(radial), OBS ^VW (westward), and OBS _VN (northward). 
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flow. To determine the sum of both effects a priori (without 
knowing the differential rotation), we decided to use the aver- 
age, (v(0 ± dx,y)) t -dx=io P x;±y, of a 20-pixel wide stripe along 
the meridian of (v(x, y)) t with respect to the distance, r, to the 
disk center as shown in Fig. [2^. Assuming radial symmetry, 
the curves of the northern and southern hemispheres are aver- 
aged. The resulting (v(0, y)) t is fitted by a 8th-order polynomial, 
v c iv(y). The curve runs from +50 m/s at the disk center to -60 m/s 
at some 700" and strong redshifts exceeding 300 m/s at the so- 
lar limb. This CLV is in good accordance with recent studies of 
the CLV of the Fe i line (e.g., |Balthasar|[T985l |Cruz Rodnguez| 
et al.|201l) . As the up- and downturns, which are caused by the 
supergranular migration at higher latitudes, are on the order of 
some 30 m/s, we neglect the impact of the meridional flowj^jon 
v c iv(y) and further calibrations. We use the polynomial v c \ Y (y) to 
create a radially symmetrical CLV model, v c i v (r), as displayed in 
Fig#- 

Since the convective blueshift is the result of the bigger 
impact of hot upwards moving granules than the cool, down- 
wards directed, intergranular fraction on a spatially averaged so- 
lar surface, the line shift depends on the observed spectral line 
and the line of sight. According to calculations based on the 
Liege atlas (Delbouille et al. 1990 ), the convective blueshift of 
the Fei spectral line at 6173.3 A in the center is Vb = -305 m/s. 
For observational studies in Sec.[3j we have to offset v c i v (r) by 
v ff ~-350m/s to obtain the convective blueshift, Vb(r), with its 
CLV. 

After the calibration of v TOt (d, /, Bo) and v res (x, y) in Sects. [23 



and |2.4| we can check the accuracy of the data calibration by 
computing the CLV velocity map a posteriori: 



v c \v(x,y) = (v(x,y)) t - v mt (d, /, B ) ■ 



s(x,y) 



Then, we compare the CLV in the meridional and 
equatorial directions, i.e. (v c iv(0±dx,y))dx=wpx;+y and 
(v c iv(X 0±dy))d y =iop\;±x- The curves are displayed in Fig. [2b 
and closely resemble each other. The velocities were offset 
by v ff to achieve the run from Vb(cos# = 1) = -305 m/s to 
Vb(cos# = 0.8 . . . 0.6) ^-400 m/s (preferred spot location) and 
lower velocities at the solar limb. While the equatorial direction 
exhibits minor fluctuations, the meridional direction displays 
stronger fluctuations owing to the residuals of supergranules. 
The difference between both curves shown in the upper part of 
Fig. |2j: reveals an rms smaller than 10 m/s in low latitudes. The 
dashed green curve in Fig. |2j: results from synthetic profiles 



for Fei 6173.3 A based on COBOLD simulations (Freytag 
et al.||2012[ |Beeck et aLl|2012| ). The synthetic curve differs by 



less than 100 m/s, so is in good qualitative agreement with our 
measurements. 



Fig. 2. — a) Average CLV, (v(0,y)) t , along the meridional cen- 
terline with an 8th-order polynomial fit, v c \ Y (y) (red curve) in 
m/s against the distance y to the disk center in arcsec. - b) 
Calibrated, radially symmetrical CLV model, v c i v (r), with axis 
given in arcsec. - c) Convective blueshift, Vb(r), in v c i v along 
the meridian (red solid curve) and equator (blue solid curve) and 
their difference (black solid curve) displayed in m/s against cos 6 
with heliocentric angle, 6. A synthesis of the spectral line from 
COBOLD simulations (green dashed curve) is added. The dotted 
black box refers to the preferred sunspot location. 



2.3. Differential rotation: v YOt (d, /, B ) 

Since the rotation velocity is differential on the solar surface and 
depends on the line of sight (LOS), we have to generate a spe- 
cific, axially symmetrical solar rotation model, v TOt (d, /, Bo), with 
respect to the distance, d, from the meridian, the latitude, /, and 
the inclination, Bo, of the rotation axis. 

We start out with the time-averaged velocity map reduced 
by the radially symmetrical component, i. e. (v(x,y)) t - v c i v (r), 
displayed in the right hand panel of Fig. [T] We determine the 
absolute latitudinal velocity, v rot (/) for Bq - 0°, by performing 



This deformation along the meridian contains the center-to-limb 

variation (CLV) of the convective blueshift and the meridional 35° ( [Hathaway & Rightmire|2010 Komm etal.|2011| ) 



3 A large-scale axisymmetrical flow directing from the equator to the 
poles on the solar surface with velocities up to 20 m/s at latitudes around 
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Fig. 3. — Absolute rotation velocities, v rot (/, #o = 0), in m/s 
(black curve) on the solar surface for latitudes up to / = 
±80° fitted by (red curve) and an equatorial maximum of 
Q(0) = 1992 m/s. 



panel), generates the LOS velocities for every position on the 
solar disk according to the annually varying Bq = ±7.27°. 
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Residual velocity map, v res (x,y), with axis given in 



linear regressions for the LOS velocities along latitudinal cuts, 
v(d/Ri), with Ri being the distance of the solar limb from the 
meridian: v rot (d, /) = v rot (/) • (d/Ri). Figure [3] shows the latitudinal 
dependence, v rot (/). 

To compare our rotation model with literature values, we fol- 
low [Stix] ([2002]) and approximate v rot (/) by 

Q(/) = A + B • sin 2 / + C • sin 4 / . (2) 

The approximation, Q(/), is plotted in Fig. [3] We obtained a 
differential decrease to higher latitudes. The equatorial maxi- 
mum, A= 14.10 ±0.037day, corresponds to D(0) = 1992 m/s. 
The coefficients B = -9.0° /day and C = -2.5°/day have stan- 
dard deviations on the order of 0.1° /day. Since the result de- 
pends on the method of measurement (see: Stix 2002 ), we com- 
pared ou r rotational veloci ty with the Doppler shift measure- 
ments of |Snodgrass| dl984| ) who obtained an equatorial rate of 
A = 14.05 ± 0.01° /day or £1(0) = 1975 m/s. This is in good accor- 
dance with our findings. However, compared to Snodgras coef- 
ficients, B- -1.49° /day and C- -2. 61° /day with deviations on 
the order of some 0.1 °/day, we obtained a stronger decrease in 
the rotational velocity to higher latitudes. 
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Fig. 4. — Computed differential rotation maps, v rot (d, /, Bq = 0°) 
(left panel) and v YOt (d, /, B^-l 2T) (right panel), with 
maximum velocities v rot (/ = 0, Bq = 0°) = 1992 m/s and 
v ro t(/ = 0,2?o = 7.27°) = 1976 m/s, spatially displayed in arc- 
sec. 

The differential rotation model, v YOt (d, /, B Q ), displayed as an 
example in Fig. g] for B = 0° (left panel) and B = 1 IT (right 



2.4. Residual velocity map: v res (x, y) 

Augural annular patterns (see Fig. [TJ indicate a systematical, 
nonsymmetrical residual, v res (x,v), that contains instrumental 
effects and other systematic flow fields of smaller magnitude. 
Subtracting v c i v and v rot from the 61 day average, (v) t , we ob- 
tain the residual nonsymmetrical velocity map: 



Vresfe y) = (v(x, y)) t - v dv (r) - v rot (d, /, B = 0°) 



(3) 



which is displayed in Fig. [5] Our result agrees well with the 
HMI c alibration studies of Howe et al. (201 1) and |Centeno et al.| 
( [201 1] ). Several large circular patterns and smaller interference 
rings were checked as spatially fixed with velocities varying be- 
tween ±150 m/s. The residual map is subtracted a priori for all 
Doppler maps in Sec. [3] Changes in time are expected to have 
small amplitude so will not affect our analysis of horizontal flow 
velocities. 



3. Data selection and analysis 

Overview: In this section, we analyze the sunspot moat flow 
and Evershed flow in the solar photosphere by the calibrated 
720s Doppler maps recorded by HMI. We have selected circu- 
lar sunspots and constructed 3h time averages (out of 15 such 
Doppler maps) of the field of view (FOV). By, doing this, were- 
duce the amplitudes of granulation and oscillation (see Sec. |3.1| ). 
We analyze the flow field in the penumbra and in the surrounding 
moat cell. 



3.1. Data selection and preprocessing: 

The moat flow is coupled to the presence of a penumbra and 
follows the direction of the penumbral filaments . Thus, we 
searched for sunspots with a fully-fledged, circular penumbra 
and found 31 applicable, single sunspots in a slight deca^in 
the time between June 2010 and January 2012. The selected 
sunspots vary in size between 9 Mm and 22 Mm and are listed 
in Table |A.1| of the appendix. A contour of the outer penumbra 



Zurich Classification: Group H 
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Fig. 6. — a) 720s intensity map and b) correspondi ng 72 0s velocity map, Vdop, of AR11084 located at 19°S,48°W, respectively 
6 = 51°, recorded 2010 July 6 at 01:36 UT (see Table A.l No. lb). An arrow is pointing to the disk center. The penumbral contour 
(white line), the theoretical, circular spot shape (black ellipse), and foreshortened circles with rising radii are displayed as well. The 
FOV size is 301 p x 150"). - c) Aligned 3h average map, V3h, of AR11084. - d-h) Masked versions of spots No. lb, 4a, 6a, 
7a, 9a (Table A.l ) at a FOV size of 201 px 100") and white ellipses with R$ + x L=5, 10,20,/^ • The velocities are displayed up to 
±1000m/s. 



for each spot is determined by means of an intensity threshold 
for a spatially averaged image. The center of a sunspot is deter- 
mined as the center-of-gravity of all points inside its penumbral 
contour. 

To analyze predominantly horizontal flows with LOS 
Doppler maps, we selected the observing position of the sunspot 
center at heliocentric angles of some 6 « 50°. There, the LOS 
component of the horizontal flow is significant, while the cali- 
bration uncertainties that increase towards the solar limb are still 
small. 

Each sunspot FOV from the size of 301 x 301 px is tracked 
by its averaged center, (x c ,y c ), for 15 successive 720s intensity 
maps (Fig. [5k), while the respective Doppler maps, v^ op (x, y), are 
calibrated By 

Vdo P (X y) = VdopU, y) - Vrot(d, /, B ) - v res (x, y) + v ff (x, y) (4) 

to obtain velocity maps, VdopC*,} 7 ), (Fig. [6b) free of all system- 
atical velocity components (described in Sec. [2]), but the surface 
convection. v ff describes the offset of the convective blueshift 
and its CLV (Fig. [2j0 to values around -400 m/s at the chosen 
observing position. 

The sunspot and flow properties are stable within three hours 
. To diminish 5-min oscillations and granulation, we averaged 
the tracked and aligned sunspot FOVs (201 x 201 px) of all the 
15 successive velocity maps and refer to this 3h average as V3h 
(Figs|6j:-h). Likewise, the received LOS velocities are adapted 
for an average heliocentric angle, 6. For better illustrating the 
MF, the 3h average is shown here with a masked umbra and 
penumbra. 

We studied the weekly evolution of the observed components 
by tracking 20 sunspots near the eastern limb to a second proper 
position (#^ 50°) after some 6 to 8 days and compared the results 
of both V3h analysis. As this was not possible for 11 of 31 spots, 



we achieve a sample of 51 velocity maps. Tracking three long- 
lasting sunspots across the fa r side of the Sun, they reappear at 
the eastern front side (Table |ATTJ No. 4 -> 6, 7 -> 9, 11 -> 13; 
Figs.[6^-h) which allows us to study the monthly evolution. 

3.2. Method of flow analysis 

We analyze flow fields of fully calibrated velocity maps, V3h, 
by applying a m ethod that determines the az imuthally averaged 
flow properties ( Schlichenmaier & Schmidt 2000). In this way, 
we assume axially symmetrical flow fields for circular sunspots 
with a fully developed penumbra. The foreshortening effects of 
a circular sunspot at a certain heliocentric angle, 6, cause the el- 
liptical shape as shown for AR11084 (Figs. [6^-d). To observe 
the spatial velocity dependency of the flow fields, we use an au- 
tomatic procedure that generates ellipses as foreshortened cir- 
cles (with radii, r, in px) based on the heliocentric angle and 
rotated according to their position on the disk. The LOS veloci- 
ties, v L0S (r, 0), along the ellipses are read out according to their 
circular angle, 0, as shown in Fig. [7] for distances r = R$ + 5 
(upper panel) and r = R$ + 10 (lower panel) from the spot cen- 
ter. The sunspot radius is automatically determined by the min- 
imal deviation of the ellipse to the penumbral contour for every 
720s intensity map and averaged for the 3h period. A manual 
fine adjustment was added for irregularities. The LOS velocities, 
v L0S (r, 0), were fitted by a sine function: 

v L0S (r, 0) = vj# s (r) • sin + v^ os (r) . (5) 

The amplitude, v^ s , is the LOS component of the horizontal 
flow velocity 

vn,o(r) = v^ s (r)/sin^ , (6) 
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vV 0S , is the vertical 



at the heliocentric angle, 6. The axial offset, 
flow component. 

3.3. Moat flow 



We analyze the moat flow in the vicinity of the s unspo ts for the 
fully calibrated velocity maps, V3h (listed in Table | A. 1 }. 



3.3.1. Moat flow velocity 

The moat flow is an outwardly directed flow, which is predom- 
inantly horizontal. We present the flow field analysis for the 
velocity map, v 3h , of AR11084 (Figs.^-d; Table [XT] No. lb). 
The sunspot in the FOV center with an average spot radius of 
Rs = 31 px « 11.2 Mm at = 51° exhibits a distinct moat flow 
surrounding the sunspot with LOS velocities up to -1000 m/s 
for the blueshifted side facing the disk center and +1000 m/s for 
the redshifted limb side indicated by an arrow. The outer vicinity 
of the sunspot features supergranular flow cells. 
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Fig. 7. — Angular LOS velocities, v L0S (r, 0), in m/s along the el- 
liptical projection for = 51° of circles with radius r = 7? s + 5px 
(left panel) and r = R$ + lOpx (right panel) fitted by a sine curve 



(green, solid line) according to Eq. |5| The amplitude, vjp q (r), 
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and offset, v^ us (r) are marked by red and blue dashed lines. 



Fig. 8. — Upper panel: Horizontal flow velocities, V||,o(r), in m/s 
for the distance, r, from the spot center in Mm. Error bars rep- 
resent the standard deviation. Dashed vertical lines mark the 
boundaries between umbra (U), penumbra (PU), and moat flow 
region (MF) and the quiet Sun (QS). — Lower panel: Horizontal 
LOS velocity, v^ s (r), in m/s of the MF (green curve) and the 
quiet Sun (red curve) against the distance from the spot in pixels 
with <x rms (see Eq.[7} as error bars. 



Comparing the velocity components of the MF in a distance 
of 5 px and lOpx from the sunspot (see Fig. [7]) 

v\° s (R s + 5px) = -792 ± 17 m/s 
v\° s (R s + lOpx) = -555 ± 12 m/s 
v^ os (tf s + 5px) = -266 ± 12 m/s 
v^ os (tf s + lOpx) = -259 ± 9 m/s , 

the most prominent result is the decrease in v|j^ s by more than 
230 m/s within these 2 Mm. The upper panel of Fig. [8] shows 
the continuous run of the horizontal flow velocities, V\\ t o(r), at 
6 = 51° (see Eq.[6]), from the spot center to a distance of 25 Mm. 
The MF has a strong monotone decrease with increasing dis- 
tance from 1150 m/s just outside the penumbra to 600 m/s after 
additional 3 Mm and 250 m/s after 6 Mm. 

The vertical flow component of MF and EF depends on the 
choice of the rest frame. Our rest frame is determined by the 
convective blueshift of nonmagnetic convection. Since the MF 
region is not completely free of magnetic fields, it may change 
the convective blueshift by an unknown amount. Therefore, no 



accurate analysis of the vertical MF component, v^ os (r), can be 
given. Also, the inclination of the MF could not be reliably deter- 
mined owing to the dependence on the vertical flow component 
( Schlichenmaier & Schmidt 2000 ). However, on the assumption 
that no effect is caused by magnetic fields, we can subtract the 
convective blueshift Vb(r)^-400m/s from v^ os (r) and obtain 
absolute vertical velocities of v^ os (r) = . . . 200 m/s. According 
to that, the slight redshift indicates a downward-directed MF. In 
comparison with horizontal MF velocities up to s = 1000 m/s, 
the vertical velocity is an order of magnitude less, and therefore 
it has a minor impact on the absolute flow velocity. 



3.3.2. Moat flow extension 

To detect the size of the MF region, we define the criteria 
v^[? s < <r rms to indicate the end of the detectable MF, whereas 
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the root mean square, <x rms , is calculated for each distance, r, 
the spot by 



to 



sir) 



1 n 2 

-•£(v^(r, fc) -/(ft)) 

i=l 



(7) 



In the lower panel of Fig. [8] we display <x rms (r) as error bars 
to the horizontal LOS velocities. When <x rms (r), i.e. the velocity 
fluctuation around the sine fit, f((f>i), with circular angles, (pi (see 
Fig. [Tj), exceeds the fitted amplitude, v^ s , we define this dis- 
tance as the end of the MF region. For the case under consider- 
ation, we obtain an MF, which extends 6.5 Mm into the vicinity 
of the sunspot. This study of the MF extension highly depends 
on cr rms , which is about 180m/s for the 3h average and which 
becomes smaller for longer periods. But, these longer periods 
cannot guarantee the stability of the observed components. We 
therefore believe that our method delivers a good estimate for 
the MF extension. 



3.4. Evershed flow 

Also the well-known penumbral Evershed flow is a predomi- 
nantly horizontal flow. In our analysis we assume it as axially 
symmetrical. In Fig. [8] (upper panel), we display the radial de- 
pendence of the horizontal flow component of AR 11084 (see 
Fig. |6j:) based on V3h. The inner penumbral boundary is deter- 
mined as the average extension of the umbral darkening in all 
720s intensity maps. We observe EF velocities, V|| ? o(r), increas- 
ing from the inner penumbra to a maximum value of 2210 m/s at 
a distance r = 8.5 Mm^0.76-7? s followed by a decrease to about 
1150 m/s at the outer penumbra. Based on our analysis method, 
some penumbral filaments will reach into the defined MF re- 
gion or end earlier. Therefore, the velocities at the transition re- 
gion have impacts on both flows, and therefore the transition is 
smooth as seen in Fig. [8] But, following the run of the curve, 
a kink between the radially decreasing EF and MF velocities is 
evident at the outer penumbral boundary . The mere existence of 
the kink is a hint at the separation of both flow regions . 



4. Results and discussion 

Overview: In this section, we compare the analysis results (see 
Sect. [3]) for all 51 velocity maps V3h listed in Table |A.l| We dis- 
cuss the interactions between the sunspots and their moat flows. 
The weekly evolution of the observed components is studied for 
a sample of 20 sunspots, whereas the monthly evolution is an- 
alyzed for three sunspots. We compare our findings with recent 
studies and analyze the relation of MF properties with sunspot 
properties. Our discussion is focused on whether the EF and MF 
are related and on what we can learn about the physical origin of 
the MF. 



4.1. Horizontal moat flow velocity: v M f(^) 



In the appendix, Figs. |A.l| and |A.2| we display the radial depen- 
dence of the horizontal flow component for all 51 sunspot maps. 
Since this is the first time that the EF and MF are analyzed co- 
herently for more than only a few spots, it is very remarkable to 
find that all our spots behave very similarly. The radial depen- 
dence is qualitatively the same for all spots. At first glance, only 
the sunspot size and the maximum of the EF velocity vary. As 
described in Sec. |3.4| the kink between the velocity decrease of 
the EF in the penumbra (PU) and the MF at the spot boundary 
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Fig. 9. — Horizontal MF velocities, vmfM, in m/s of 51 ob- 
served 3h velocity maps against the distance from the spot cen- 
ter in spot ratio (upper panel) and the distance to the sunspot in 
Mm (lower panel). The calculated average (black bold curve) is 
drawn. 



is recognizable for all s unsp ots in the radial dependences dis- 
played in Figs. A.l and A.2 In the next step, we examine the 
flow properties of all spots in detail and compare them. 

The moat flow was identified for all selected sunspots as 
a predominantly horizontal, axially symmetrical flow starting 
just outside of the penumbra with velocities depending on the 
distance, r, to the sunspot. Figure [9] displays the strong de- 
cay of the MF velocities, vmf(^) = V||,o(^) for all 51 analyzed 
MF regions based on V3h. The upper panel shows the decrease 
with respect to the distance from the spot center in the ra- 
dius ratio, r = ^j^. The lower panel plots the MF velocities 



against the distance, f = r = 0, . . . , 14 Mm, from the penum- 
bral boundary. The monotone decrease from maximum veloc- 
ities of vmf(7* = 0) = 800 . . . 1200 m/s just outside the sunspot is 
largely similar for all sunspots. The average run, <VMF(^))5ispots> 
starts with 1000 m/s just beyond the penumbral boundary at 
r = Mm and decreases to approx. 500 m/s within r = 4 Mm or 
r= 1.3 -7?s- The average velocity falls below the typical rms val- 
ues, <r rms « 1 80m/s, at r = 9.2 Mm or r = 1 .7 R$ . Consequently, 
the MF velocity highly depends on the location within the MF 
region. We compare both panels of Fig. [9] and notice that the 
runs according to the radius ratio in the left hand panel are 
spread wider than according to the sheer distance from the spot. 
Consequently, the MF appears to be independent of the spot size. 
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This result is in line with recent studies describing the MF 
velocity on the surface ranging betwe en vmf = 500 . . . 1000 m/s 
(Balthasa r et al.|19 96; Rimmele 1 1 997 1 [Vargas Dommg uez et al. 



2007), whereas the continuous velocity decrease with increasing 
distance from the spot boundary has not been measured before. 



4.2. Maximum Evershed flow velocity: v EF 

The maximum EF velocity, v EF = m<2x [v^o (R\j < r < Rs)], 
for all observed penumbrae is in the range of 
v EF = 1830 . . . 3000 m/s (Fig. [Kj) with an average of 2325 m/s. 
Since we measure azimuthally averaged flow speeds of 3h 
averages, our maximum EF velocities are slightly lower 
than the maximu m velocities of up to 3...4km/s found 
in recent s tudies ([Sc hlichen maier & Schmidt] |2000[ [S hine 



4.3.3. Correlation of spot size and EF velocity 

To yield the correlation of the sunspot properties, i. e. the max- 
imum horizontal EF velocity between vef = 1830 . . . 3000 m/s 
and the sunspot radius, we plot both components against each 
other in Fig. 10 1 and fit a linear regression by considering all 



values. We obtain a slightly positive trend for bigger sunspots 
harboring higher maximum EF velocities. Following the regres- 
sion, sunspots with R$ = 10 Mm would have maximum EF ve- 
locities of 2200 m/s whereas sunspots with the double size yield 
maximum EF velocities that are almost 500 m/s faster. Although 
there are some outliers between 16-17 Mm, even the correlation 
coefficient, C(<2mf, ^s) = 0-46, suggests a positive link of the EF 
velocity with the spot size. 



et al. 



1994| |Rouppe van der Voort||2002[ |Bellot Rubio et al.| 4.3.4. Correlation of Evershed flow and moat flow 



2003a). The location, tef, of the extreme EF velocity can by 
found at 0.65 Rs < tef < 0.87 Rs, with an average ratio of 
r E F = 0.78-7? s which is not as far in the outer spot boundary 
as r E F = 0.8...0.9-fl s flTritschler et~aLl [20041 |Franz| |20Tfl ). 
The velocity run within the penumbra is largely similar for 
all sunspots with some differences in the skewness as seen in 
Figs.[Al]and[Al2 



4.3. Correlations 

In the following section, we analyze the correlations between 
v M f, <2mf, vef, and Rs based on the 3h average sample. 



4.3.1 . Correlation of spot size and MF velocity 

In order to work out the interaction of the sunspots with their 
moat flows, we plot the MF velocities, VMF = (vMF(^))r=o,...,5px, 
against the sunspot size, Rs, as it can be seen in Fig. 10 \ with a 
nonsignificant rise of the linear regression (red line) and a cor- 
relation coefficient of C(vmf, ^s) = 0-28. Hence, we find that the 
MF velocity is in average around 1000 m/s and does not depend 
on the spot size. Since we measure close to the outer sunspot 
boundary the indicated slight increase may be due to a stronger 
EF in bigger sunspots (see Sec.|4.3.3]and|4.3.4"}. 



4.3.2. Correlation of spot size and MF extension 

According to Eq. [7] the MF extensions, <2mf, were determined 
for all sunspots with a boundary criteria of <x rms ^180m/s 
and range between <2 M f = 5 . . . 15 Mm for the spot sizes of 
Rs = 8.6 ... 21.2 Mm, cf. Fig. 10 3. Since the linear regression 
and the correlation coefficient of C(<2mf,^s) = 0.06 indicate 
uncorrected components, we average the MF extensions to 
(^mf )s\ spots - 9.2 Mm. Comparing our results to other studies 
(e.g.,[Brickhouse & Labonte 1988 ) the MF extension lies below 
the mentioned #mf = 10 ... 20 Mm for small sunspots and far be- 
low the double spot radius for larger spots. In our studies, the MF 
extension is independent of the spot size and therefore is in the 
order of the sunspot size only for small but not for bigger spots. 
This finding goes in line with studies by Sobotka & Roudier 
(2007]). In summary it can be stated that we find four indica- 
tions that the MF is not correlated to the spot size: (i) the kink 
of V|| 5 o (r) at the outer spot boundary (e.g. Fig. [8j upper panel), 
(ii) the smaller spread of all 51 MFs by merely aisplaying them 
with their radial distance from the spot in Fig. [9] (iii) the trend 
in Fig.[T0fc and (iv) the trend in Fig.[T0fr. 



As we pointed out in Sects. 4.3.1 and 4.3.2 the moat flow seems 
to be independent of the sunspot radius. In the following, we ex- 
amine the largely unknown link between the penumbral EF and 
the adjacent MF. It is therefore important to find out whether 
higher EF velocities also cause higher MF velocities or a bigger 
MF extension, which would argue for a partial drive of the MF 
by the EF. Due to Figs.[T()]:-d which plot the MF velocities, vmf, 



and MF extensions, <2mf, against the EF velocities, v E f, with cor- 
relation coefficients of C(vmf, vef) = 0.07 and C(<2mf, v ef) = 0.01, 
we draw the conclusion that the EF has no impact on the MF. 
This means that there is a greater decrease in speed for spots 
with higher maximum EF velocities toward the outer penumbral 
boundary (see Figs. A.l and A.2 ). Although the existence of the 
MF is directly coupled to the existence of a penumbra and its EF, 
the observed properties indicate the MF and the EF are indepen- 
dent sunspot flows. 

4.4. Weekly evolution 

To study the temporal evolution of the analyzed components, 
i. e. Rs, v mf, ^mf, and vef, we have tracked 20 sunspots (see 
Table A. 1 ) from the eastern to the western limb sides and per- 
formed the analysis of V3h for a second time at 6 « 50° . According 
to the latitude and the differential rotations, the observing dates 
range between 6 to 8 days. 

The weekly changes are in the range of AR S = -3 ... 1 Mm, 
Av E F = ±500m/s, Av M F = ±200m/s and Aa M F = -8 . . . 3 Mm. 
Figures [lOf-g plot the changes, Avmf and A^mf, of the MF ve- 
locities and extensions against the change of the spot radius, 
ARs, with linear regressions according to the correlation co- 
efficients, C(Avmf,A7?s) = -0.41 and C(A<3 M f, AR s ) = -0.23, 
indicating slightly negative trends. Following this, the MF ve- 
locity increases for strongly decaying sunspots and tends to 
decrease slightly for small ARs. The MF extension, #mf, de- 
creases for the majority of the sunspots with some exceptions 
for strongly decaying sunspots, but yields no significant corre- 
lation with the sunspot evolution. Figure 10 1 plots the change 
in the MF velocities, Avmf, against the changing MF extension, 
A<2mf, with no significant correlation, C(Avmf, A<2 M f) = 0.10. 
Because they are C(Av E f, A7? s ) = -0.06, C(Av M f, Av E f) = 0.23, 
C(A<2mf, Av E f) = -0.18, all other correlations are insignificant 
and yield no further impact. According to these results, we 
would draw the conclusions, that a strong sunspot decay leads 
to an additional drive of the moat flow by accelerating its veloc- 
ity and sporadically expanding its outreach, whereas it has no 
impact on the EF velocity. 
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Fig. 10. — Correlations of sunspot and MF properties for all 51 sunspot maps (a-e), V3h, and their evolution within one week for 
20 sunspots (f-h) indicated by linear regressions (red lines) and the coefficients, C(-, ■)■ 

- a) MF velocities, v M f, in m/s and b) MF extensions, <2mf, in Mm against the sunspot radii, 7? s , in Mm with the average values, 
(vmf)si = 982 m/s and (<2mf>51 = 9.2 Mm. - c) MF velocities, vmf, in m/s and d) MF extensions, <2mf, in Mm against the maximum 
EF velocities, vef, in m/s. - e) Maximum EF velocities, vef, in m/s against Rs in Mm with the average value, (vef)si = 2324 m/s. - 
f) Weekly changes, Avmf, in m/s and g) Aa M F in Mm of the moat flow against AR S in Mm. - h) Weekly changes of MF properties, 
Avmf, in m/s against A^mf in Mm. 



4.5. Monthly evolution 



5. Conclusion 



Tracking three long-lasting sunspots across the far s ide o f the 
Sun, they reappear at the eastern front side (Table |A.1| No. 
4 — > 6, 7 — > 9, 1 1 — > 13; Figs. |6^-h). This allows us to study the 
long-term evolution for one month in order to crosscheck the re- 
sults of the weekly evolution. The changes in sunspot size, ARs, 
the MF velocity, Avmf, the MF extension, AaMF, and the max- 
imum EF velocity, Avef, as listed in Table T] are based on the 
mean difference between the first and second appearances. 



Table 1. Monthly evolution for three long-lasting sunspots. 



5.1. Conclusions on the moat flow 



Spot-No. 


AR S 


Aa MF 


Av MF 


Av EF 


(TableCO) 


(Mm) 


(Mm) 


(m/s) 


(m/s) 


4 J 6 f 


-2.8 


-3.1 


+77 


-65 


7^9 


-6.1 


+6.6 


+ 173 


-455 


11 13 


-5.2 


-5.4 


+64 


-115 



Coupled to the strong decay in sunspot sizes, the MF veloc- 
ity significantly increases, whereas the EF velocity shows a de- 
crease of several hundred m/s. These trends are in line with the 
correlations for the weekly evolution . The size of the MF region 
shows no unique trend, but the huge widening AaMF = +6.6 Mm 
of Spot-No. 7 (in Table[T]) by more than 6 Mm could underline an 
additional outflow of plasma over the moat cell due to a strong 
sunspot decay. 



Mey er et aL] ( |1979| ) suggested that a sunspot is embedded in a 
supergranular cell. If one puts MF cells in a context with super- 
granules, one should realize that they have a diameter between 
30 and 60 Mm, i.e., are up to twice as large as supergranules. 
Also, the MF velocities are more than twice as high as super- 
granular velocities (e.g., Brickhou se & Labonte|1988| ). As a re- 
sult, the MF needs a driving mechanism that is distinct from the 
driving mechanism of normal supergranules. Nye et aL| ( 1988) 
modeled the MF being driven by a surplus gas pressure beneath 
the penumbra, arguing that the lack of radiative cooling beneath 
the penumbra generates a surplus of heat, hence gas pressure. 
In these models the MF velocity depends on the depth of the 
penumbra, but not on the size of a spot. At that time this posed 
a problem, as it was commonly believed that the MF extension 
scales with the spot radius. Later, Sobot ka~& Roudier] ( |2007| ) 
provided evidence that there is no correlation between the MF 
and the spot radius. Our investigation based on Doppler shift 
measurements rather than on local correlation tracking, indepen- 
dently confirms these later findings. The evidence therefore sup- 
ports the moat flow model in which the driving forces are due to 
surplus gas pressure beneath sunspots. 

A moat flow that is driven from beneath the penumbra would 
naturally push away the granules from the spot as seen with local 
correlation techniques. This implies that the MF is present in the 
deep photosphere. In contrast, the Evershed flow or the fraction 
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of it that extends from the penumbra outwards into the moat is 
present in the magnetic canopy that surrounds the sunspot in the 
mid and upper photosphere (Rez aei et al.|20 06). In the immedi- 
ate surroundings of the sunspot two types of flows exist: (1) the 
(magnetic) EF that partially continues in the magnetic canopy, 
which ascends outwards from the mid-photosphere at the spot 
boundary; (2) the (largely nonmagnetic) MF in the deep photo- 
sphere and beneath. 

We consider the moving magnetic features (MMFs) to be 
distinct from the moat flow. They migrate into or are advected by 
the MF. MMFs are associated with inclined magnetic field lines 
(relative to horizontal) that reach up into the higher atmospheric 
layers. These MMF field lines possibly originate in the sunspot 
and either witness the decay of sunspots (uni-polar MMFs) or 
are due to some waves that propagate outwards (bipolar MMFs), 
(see e.g., |Sainz Dalda & Martinez Pillet||2005[ |Sainz Dalda &| 
Bellot Rubio 2008a b; Schlichenmaier 2002 ). Because of their 
vertical configuration, outwardly migrating MMFs can also be 
detected in hi gher atmospheric layers as reported by Sobotka & 
|Roudier| ( [2007l ). 

5.2. Summary 

We calibrated HMI velocity maps such that systematic absolute 
errors are below 150 m/s on an absolute scale . The synoptic 
CLV of the convective blueshift was measured in good accor- 
dance with previous findings and the theoretical synthesis curve. 
We find a maximum velocity of the solar rotation of 1992 m/s, 
which agrees with previous measurements and reveals no signif- 
icant impact by stray light. The instrument al artifacts were iden - 
tified in line with HMI calibration studies (Centeno et al.|20lT] ). 
These artifacts make the meridional flow amplitude too small to 
be measured. 

We analyzed 3h time averages of sunspot flow velocities. We 
constructed 51 velocity maps of 31 sunspots. The flow in and 
around circular sunspots with a fully developed penumbra was 
analyzed using azimuthal averages, thereby assuming axial sym- 
metry. In both, the MF and EF, the horizontal velocity compo- 
nent was dominant. The vertical flow components of MF and EF 
were determined to be small. Since the exact amount of convec- 
tive blueshift is unknown, we do not address the sign of vertical 
flows. 

The radial dependence of the velocity fields was observed to 
be similar for all 5 1 sunspot maps. The analysis of the EF yielded 
results that are consistent with recent studies. A higher EF veloc- 
ity was detected for bigger sunspots, and the EF velocity scales 
with the spot size. The MF is a convective outflow with radially 
decreasing velocity. The MF velocity turned out to be similar in 
size and independently of the spot size. Also, the MF extension 
is not correlated with any sunspot property. As MF and EF prop- 
erties turn out to be uncorrelated, we inferred that the driving 
mechanisms of the two flows should be distinct. Therefore, we 
favor the MF model of |Nye et alT] ( [T988j ) : The MF is driven by 
surplus gas pressure beneath the penumbra. 

We find a tendency toward increasing MF velocity for 
strongly decaying sunspots on time scales of one week (spot ro- 
tates from the east to west limbs) and one month (spot reappears 
on the east limb). It is beyond the scope of this paper to put this 
into the context of a model for decaying sunspots. 
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Appendix A: Analyzed sunspots 

Table A.l. List of the 31 numbered, observed sunspots. The corresponding NOAA number of the active region (AR) is given. 
The number, date, and starting time of 3h observation, as well as the average spot location, the spot radius, Rs, (in Mm), and the 
heliocentric angle, 6 (in °), are listed. The spots with No. 4+6, 7+9, and 11 + 13 describe the same recurrent sunspots, respectively. 



No. 


AR 


Obs.-No. 


Date 


Time 


Location 


Rs [Mm] 


6 [°] 


1 


11084 


a 


2010-06-29 


00:00 


S19E45 


11.5 


49 






b 


2010-07-06 


00:00 


S19 W47 


11.2 


51 


2 


11092 


a 


2010-08-06 


20:00 


N13 W44 


15.1 


45 


3 


11101 


a 


2010-08-27 


00:00 


N13 E45 


11.5 


46 






b 


2010-09-03 


08:00 


N12 W48 


11.9 


49 


4(->6) 


11108 


a 


2010-09-25 


00:00 


S30 W36 


17.3 


49 






b 


2010-09-25 


10:00 


S30 W40 


17.3 


52 


5 


11113 


a 


2010-10-16 


00:00 


N17 E48 


12.6 


49 






b 


2010-10-23 


10:00 


N16 W48 


10.4 


48 


6(<-4) 


11115 


a 


2010-10-17 


18:00 


S28 E40 


15.1 


52 






b 


2010-10-24 


04:00 


S29 W40 


13.7 


52 


7(->9) 


11131 


a 


2010-12-05 


00:00 


N30 E39 


21.2 


49 






b 


2010-12-11 


08:00 


N30 W38 


19.8 


47 


8 


11133 


a 


2010-12-14 


04:00 


N14 W49 


10.1 


51 


9(<-7) 


11140 


a 


2011-01-02 


20:00 


N32 E40 


15.1 


51 






b 


2011-01-08 


20:00 


N33 W35 


13.3 


47 


10 


11147 


a 


2011-01-18 


16:00 


N25 E41 


12.6 


45 






b 


2011-01-24 


12:00 


N24 W43 


11.9 


49 


11 13) 


11195 


a 


2011-04-28 


10:00 


S17 W44 


15.1 


48 


12 


11203 


a 


2011-05-01 


12:00 


N18 E45 


12.2 


47 






b 


2011-05-08 


09:00 


N17 W47 


10.8 


50 


13 11) 


11216 


a 


2011-05-18 


12:00 


S15E47 


10.4 


48 






b 


2011-05-25 


14:00 


S15 W46 


9.0 


47 


14 


11251 


a 


2011-07-13 


15:00 


N16 E49 


12.2 


49 






b 


2011-07-20 


20:00 


N17 W46 


9.4 


46 


15 


11260 


a 


2011-08-02 


18:00 


N19 W46 


17.6 


50 


16 


11277 


a 


2011-08-27 


00:00 


N17 E47 


9.0 


51 






b 


2011-09-03 


09:00 


N19 W46 


9.7 


50 


17 


11287 


a 


2011-09-06 


00:00 


S31 E37 


12.2 


49 






b 


2011-09-11 


18:00 


S29 W39 


9.0 


52 


18 


11305 


a 


2011-10-04 


20:00 


Nil W48 


14.8 


50 


19 


11312 


a 


2011-10-07 


00:00 


N22 E48 


17.6 


49 






b 


2011-10-14 


03:00 


N23 W43 


17.3 


45 


20 


11314 


a 


2011-10-19 


00:00 


N26 W42 


16.6 


46 


21 


11317 


a 


2011-10-14 


00:00 


S27 E42 


9.4 


51 






b 


2011-10-20 


06:00 


S27 W41 


9.4 


51 


22 


11338 


a 


2011-11-10 


17:00 


S12 W48 


12.6 


53 


23 


11340 


a 


2011-11-07 


00:00 


S09 E49 


10.1 


50 






b 


2011-11-14 


10:00 


S08 W49 


9.4 


50 


24 


11342 


a 


2011-11-15 


00:00 


N17 W47 


12.6 


48 


25 


11343 


a 


2011-11-16 


00:12 


N28 W39 


8.6 


44 


26 


11355 


a 


2011-11-27 


21:00 


N14 W45 


9.7 


47 


27 


11356 


a 


2011-11-27 


21:00 


N16 W23 


10.8 


30 


28 


11366 


a 


2011-12-05 


00:00 


N17E51 


9.7 


53 






a 


2011-12-12 


08:00 


N18 W47 


9.0 


48 


29 


11384 


a 


2011-12-23 


00:00 


N12 E36 


19.8 


37 






b 


2011-12-26 


00:00 


N13 W01 


18.7 


14 


30 


11388 


a 


2011-12-30 


00:00 


S23 E42 


10.8 


45 






b 


2012-01-05 


00:00 


S24 W35 


9.7 


40 


31 


11389 


a 


2011-12-31 


00:00 


S23 E43 


15.8 


45 






b 


2012-01-07 


00:00 


S21 W45 


16.9 


50 
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Fig. A.l. Horizontal flow velocities, V\io(r), in km/s from the sunspot center (r = OMm) to a distance of 36Mm in the quiet Sun(QS) 
for spots No. 1-17 (see Table A.l }. The boundaries of the umbra (U), the penumbra (PU), and the end of the moat flow region (MF) 
are indicated as vertical dashed lines. 
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Fig. A.2. Horizontal flow velocities, V\\ t o(r), in km/s from the sunspot center (r = OMm) to a distance of 36Mm in the quiet Sun(QS) 
for spots No. 18-31 (see Table A.l }. The boundaries of the umbra (U), the penumbra (PU), and the end of the moat flow region 
(MF) are indicated as vertical dashed lines. 
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